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Introduction 

A. Why should hydromaqneti c waves be measured in space ? 

HM waves are produced by plasma instabilities or currents 
far out in the magnetosphere and propagate to the earth's sur- 
face. For many hm waves there is reason to believe that propa- 
gation is f iel d-dlii gned , or ducted. Thus, the waves seen at the 
earth's surface provide information about plasma processes going 
on at known points in space. The wave spectrum, polarization, 
and direction of propagation can be used to infer the nature of 
these processes. For example, the existence of Pc 1 waves 
(left elliptically polarized transverse waves of 5 sec period) 
has been used to infer the occurrence of the ion cyclotron in- 
stability of energetic protons in the magnetosphere. Furthermore 
the fact that one type of Pc 1 (pearl pulsations) shows disperion 
in its dynamic spectrum has been used to determine the L shell 
on which the waves are propagating, the equatorial density of cold 
plasma on this line, and the parallel energy of the protons re- 
sponsible for the instability. Another type of Pc 1, IPDP, is 
beginning to provide information about the point and times of 
proton injection during substorms as well as details about the 
drift of these particles. Other possibilities include the high 
latitude Pc 4,5 seen near the polar cusp, or Pc 3,4 seen near 
the plasmapause. These waves may eventually be used to provide 
information about the location of these boundaries as well as 
processes occurring on them. 

In summary, hm waves are diagnostic of plasma processes in 
the magnetosphere. The location at which these waves occur on 
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the earth's surface reflects locations of important regions or 
boundaries in space. The properties of the waves reveal infor- 
mation about the plasma properties in these regions. The waves 
therefore provide the possibility of inexpensive monitoring of 
these magnetospheri c processes with ground observations. 

B. Current Feasibility of Monitoring Magnetospheric Processes 
with hm Waves 

The general goal outlined above is not presently attainable 
for a numer of reasons. Mainly this is a consequence of our lack 
of understanding of the processes of wave generation and propa- 
gation to the earth's surface. This is, in part, the major justi- 
fication for an hm wave sensor on Space Shuttle. We presently 
have almost no understanding of the effects of the ionosphere on 
the transfer of waves. For some waves the ionosphere may act as 
a resonant cavity, increasing the amplitude of the incident waves 
much above their initial amplitude. Alternatively, the ionosphere 
may duct energy away from the point of entry. Almost certainly 
the polarization of the incident waves is greatly altered by the 
ionosphere. 

These facts make it impossible to use current ground obser- 
vations to determine the point of entry of hm waves from the 
magnetosphere, or their polarization as it exists within the 
magnetosphere. Without this information, no clear inferences 
can be made about the locations or types of magnetospheric genera- 
tion mechanisms. 

A major goal of an hm wave sensor on Space Shuttle would be 
to determine the effects of the ionosphere on the transmission of 
hm waves to the ground. 



C. Possible Uses of Space Shuttle In Calibration of Ground 
Observations of hm Waves 

1 . Latitude of hm Entry 

The most straightforward use of Space Shuttle in hm wave 
studies would be the determination of the latitude of entry of 
the waves from space. To accomplish this, the Shuttle would be 
placed in polar orbit well above the ionosphere. Since the 
Shuttle orbit would remain relatively fixed in inertial space, 
the observations would be made at fixed local time. Consequently, 
it would be necessary to choose this local time on the basis of 
information about the most probable local time of occurrence of 
the phenomena to be studied. Continuous monitoring of the output 
of the sensors would provide a body of data which should include 
examples of the phenomena to be studied. 

Most hm wave phenomena are sporadic in their occurrence as 
well as localized in space. Consequently analysis of the Shuttle 
data cannot be carried out without simultaneous ground data moni- 
toring the temporal changes in hm wave amplitude on the ground. 

If hm waves are present on the ground at a given local time and 
latitude but not at the Shuttle when it passes over the ground 
station, one would conclude the waves are propagating in the 
ionosphere to the ground station. On the other hand, if the waves 
are also seen at the satellite, then the normalized wave amplitude 
versus latitude would indicate the latitude of entry of the waves 
into the ionosphere. 

In addition to normalized wave amplitude, the direction of 
wave propagation and polarization at the Shuttle might reveal 
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the latitude of hm wave entry. For example, it has been reported 
that Pc 4 waves change their polarization across a demarcation 
line. Some recent theoretical work suggests there is ample rea- 
son for expecting such a result. 

2. Transfer Function of the Ionosphere 

In some cases hm waves will be observed both on the ground 
and simultaneously in space. In such cases the incident wave 
amplitude and polarization can be compared to that observed on 
the ground. This comparison would define the transfer function 
of the ionosphere above the station. 

3. Hm Wave Propagation in the Ionospheric Cavity 

It is well known that some hm waves propagate long distances 
in the ionosphere. Pc 1 waves have been observed simultaneously 
at conjugate points of an auroral zone magnetic field line as 
well as the equator. Dynamic spectra at the conjugate points 
contain periodic structures such that when one spectrum is dis- 
placed by half the period of the structures, the dynamic spectra 
appear identical. In contrast, the equatorial spectrum appeared 
to be the sum of the conjugate spectra. Using these kind of data 
the group delay between stations can be calculated as a function 
of latitude. The results obtained correspond roughly to the 
Alfven wave velocity in the F region in the ionosphere. 

This ionospheric propagation of hydromagnetic waves has 
been studied theoretically as well. Among the question of in- 
terest are how the waves are injected into the cavity, how large 
are wave amplitudes in the cavity, how does wave polarization 
change within the duct, how rapidly are waves attenuated in the 
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Instrumentati on 

The construction of a hydromagneti c wave sensor for Space 
Shuttle presents a number of difficult technical problems. These 
problems are centered around the necessity of measuring very small 
fluctuations in the presence of an extremely large, but variable 
background field. If we accept 10 my as the desired resolution 
and ±50,000 y as the necessary dynamic range for the earth's 
field, we require a resolution of 10 x 10 /lO or one part in 10 . 
If we express this instead in binary form we have roughly: 

(l/64)/131 ,072. = l/(2®)(2^^) = 1/2^^ = 1/8.39 x 10® 

For comparison, binary analog-to-digital converters of 16-bit 

accuracy have only become commercially available in the last few 

years. As yet these devices have not been qualified for space 

applications by NASA. The problem is comparable even when we 

give up the requirement of absolute field measurement and require 

only the time derivative of the field. As we show in another 

section, the maximum rate of change in any component of the earth's 

file is 50 Try per second. For a 10 my wave amplitude the time 

_2 

derivative is 2TrfA = 27 t x 10 f. The lower limit for hm waves 

- 3 5 

is roughly 10 Hz, so we require a resolution of 2tt x 10“ /50it = 

1/2.5 X 10®. 

Actually, the foregoing argument exaggerates the problem be- 
cause we have not considered certain features of the natural geo- 
magnetic spectrum and of electronic noise in instruments. In 
particular, both electronic noise and amplitude of hm waves de- 
crease as frequency Increases. For example, typical wave amplitudes 
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_ 3 

at 10 Hz are 100 y. At higher frequencies they decrease roughly 
in inverse proportion to frequency, being about 100 my at 1 Hz. 

If we require a resolution of one part in 100 of the typical wave 

_ 3 

amplitude, then we need only 1 y at 10 Hz, so that our instru- 
ment resolution must be: 

27r(l )(10'^)/50 TT = 1/(2. 5 x 10^) 

Note that 1/2^^ = 1/(1.6384 x 10^). Thus, a 14-bit digitizer 
would be nearly adequate. As we will discuss below, there are 
similar strategies in the design of absolute instruments that 

make comparable reductions in the required resolution. 

14 15 

We emphasize that a resolution of one part in 2 -2 is 

not a simple matter to achieve in a spacecraft environment. A 

typical instrument using modern integrated circuits will have a 

.2 

maximum voltage output of ±10.24 v. This is 2 • 1024 x 10 v or 
2 X (10 v). As a rule of thumb, most spacecraft ground lines 

_ 3 

will have about 5 mv of noise in them. This is 5x10 = 0.5 

-2 -2 

X 10 - 1/2 X 10 V. Thus, expressed as a fraction of full 

scale voltage, the noise is (l/2_x 10“^)/(2-^ x 10“^) = 1/2^^. 

It would appear pointless to attempt measurement of such a signal 
15 

to 1/2 . We note, however, if the noise on ground lines could 
be reduced within the instrument to 1.25 x 10“ v = 0.125 x 10“ = 
1/8 X 10“^, we would have (l/23)/(2^^) = 1/2^^, which is roughly 
adequate as shown above. 

In the previous paragraph we have continued to ignore the 
spectral characteristics of the hm signal, spacecraft noise and the 
process of digitization. To proceed further in our discussion of 
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the feasibility of an hm wave sensor for Space Shuttle we must 
consider these. It is a well-known fact that digitization of data 
introduces white noise. The spectral density of this noise is 
given by 

Pjj = (D)2/12B (v^)/Hz 

where D is the least significant bit of the digitization in 
volts and B is the bandwidth of the quantized signal in Hz. Since 
we are considering sampled data, it can be assumed that the band- 
width has been limited to the Nyquist frequency, fj^ = l/2At, where 
At is the time between samples. 

Noise spectra of electronic instruments are generally pro“ 
portional to 1/frequency. Furthermore, quoted noise values are 
typically determined at low frequencies by recording a few minutes 
of data on chart paper. Thus, if 5 mv is a correct figure for the 
bandwidths of interest in ULF wave measurements, it implies that 
the rms power in a 1/f spectrum between 0.01 and 1 Hz is 5 mv. 

But ^ 

(rms(fL - fy))^ = /'^ 0/f) df = P^*ln (f^/fL) 

Thus 

* (rms)^/ln 

For the above assumptions 

P^ = (5 X 10"^)^/4. 60517 V.5 x 10‘® (v)^/Hz 

Thus* P^, (f) = 5 X 10“®/f (v)^/Hz. 

iAverages of many power spectra for different types of hm 
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waves fall off as 1 / (frequency ) , i.e., wave amplitude propor- 
tional to 1/f. Thus 

(rms )2 = /'' (Pg/f^) df = P 5 f 

f^ i^^l. 

or 

(rms)^ = P 3 (l/f^^ - l/f^ 2 ) 

Pg = (rms)^/(l/fL^ - 1/f/) 

Most quasi-sinusoidal wave events have band widths comparable to 

center frequency. Thus, if we take 10 y as a typical rms ampli- 
_ 3 

tude at 10 Hz, we have 

Pg = (10)^/(l/5 X 10’^)^ - 1/(1. 5 X 10’^)^) 

Thus 

Pj = 2.81 X 10 ‘® 

and 

Pg (f) = 2.8 X 10“®/f^ 

Note that this corresponds to a 10 my rms amplitude at 1 Hz. 
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A. Search Coil 

If we are considering a search coil magnetometer, we need 
the spectrum of the derivative. If B = sin wt, the 3B/9t = 

o)B^ cos cot = 2irf B cos cot. Since the 90° phase shift is not 
important in the power spectrum, we see the signal amplitude 
differs from that of the original by the factor (2irf). Thus, 
the spectrum of the derivative is given by (2-rTf) times the 
spectrum of the original signal. Actually, we must scale the 
output of the search coil instrument by a constant (k) so that 
the largest expected signal remains on scale. Thus, the natural 
spectrum must be multiplied by the constant (2irfk) . Threfore 

Pg (f) = (47T^f^k^)(P^^) = 39.478 k^ P^ 
or 

Pg (f) = 40 k^ P^ 

The constant (k) is set by the requirement that the maximum sig- 
nal SOtt y/sec does not exceed 10.24 v. Thus, 10.24 v = k (50"Tr 
y/sec) or k = 6.5189 x 10 v7(y/sec). Thus 

V (volts) = k (2Trf Bq) 
or 

V (volts) = 2 tt(6.5189 x 10"^) fB^ 

V (volts) = (0.41) fB^ 

For comparison, a typical ground search coil system for the auroral 
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zone is scaled so that the constant multiplying fB^ is 1.0. At 
mid-latitudes it is typically 10.0. Thus, the satellite system 
must be 2.4 to 24 times less sensitive than ground systems to 
avoid saturation by motion through the dipole field. 

Returning to our "typical" search coil spectrum, we have 

Ps(f) = 4Tr^k^Pg = 4 tt^ (4.1 x 10"^^ (2.81x10"^) 
or 

P^(f) = 1.86 X 10"^ 


If we now plot our typical signal and noise spectra on log 
log plots, we have the result shown below. 



FRtQCHZ.) 
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By our preceding arguments the search coil spectrum for natural 
signals is flat. In contrast, the instrument noise falls off as 
1/f. Using the typical values discussed above, we see that we would 
not be able to observe natural signals below 0.03 Hz. We note, 
however, that small changes in the typical noise or the scaling 
of the search coil change this considerably. For example, if the 
instrument noise in the band 0.01 to 1 Hz is actually 1 mv rms we 

_ 3 

would be able to observe all the way to 10 Hz. Alternatively, 

we could allow one of the sensors to occasionally saturate and 

scale the search coil such that 2-nk = 1.0. Then, even for the con- 

_2 

servative noise estimate of 5 mv, we could observe to 10 Hz. 

As we will show in a later section, there are several reasons 
why the instrument should not attempt to measure hm waves much be- 
low 0.03 Hz. Consequently, we feel that there is no need to allow 
the instrument to saturate and, furthermore, only small improvements 
in instrument noise would be required. 

Finally, let us decide what the least significant bit (D) 
should be in our analog-to-digital conversion. This, in turn, 
depends on desired bandwidth. For discussion purposes, we suppose 
that we wish to measure up to f^ = 5 Hz, i.e., sample 10 times a 
second. Further, let us require that the digitization noise power 
density be 3.3 db below the instrument noise at this upper fre- 
quency limit. Then 

Pp = (2/3) P„ (f,) = D^/12 (f^) 
or 

0 • (8P„(f,)-f,)’‘'^ ■ (8.P^)’/2 2.8 /TJ 
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Substituting, we have 

D = (8 • (5 X 10"®/fy).fy)^''^ = 6 tnv 

Thus, D should be comparable to our assumed noise over two decades 
of frequency. Thus, our resolution should be about 

5 mv/ 2.10,240 mv = 1/4046 = 1/2^^ 

If instrument noise is actually 1 mv rms, then 

D = 2.8 /2 X lO-T" ■ -vl.25 mv 

and the resolution should be about 

1.25/20,480 = 1/2^^ 

As we mentioned earlier, 12-bit digitization is trivial with 
flight qualified hardware, and 14-bit should be readily obtain- 
able at some increase in cost. 

Let us now summarize our discussion of the proposed search 
coil magnetometer. We assumed that the maximum output voltage of 
the instrument was ±10.24 volts. Next, we estimated the spectrum 
of the instrument noise to be a 1/f spectrum and 5 mv rms in the 
two decades 0.01 to 1 Hz. We then scaled the Instrument such that 
the maximum expected rate of change of field (SOtt y/sec) would 
just saturate the instrument. We chose the least significant bit 
such that digitization noise power d;ensity was 2/3 of the instru- 
ment noise power density at 5 Hz. With these assumptions, we found 

1 2 

that the allowable resolution was only one part in 2 . We found 

further that if rms instrument noise was reduced by a factor of 5, 
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the allowable resolution is correspondingly increased to 1 part 
1 n z 

As a crude estimate of the desired lower limit for the 

2 

natural hm spectrum, we chose a 1/f spectrum corresponding to a 
10 my wave amplitude of 1 Hz bandwidth centered at 1 Hz. These 
assumptions led to a flat spectrum output from the search coil 
magnetometer. Because the noise spectrum rises toward lower fre- 
quencies, it is impossible to maintain constant ratio of signal 
to noise as a function of frequency. In fact, scaling the in- 
strument to avoid saturation due to motion in the dipole field 
and accepting 5 mv rms noise implies instrument noise exceeds our 
desired lower limit of the hm wave spectrum below 0.03 Hz. By 
reducing the rms instrument noise to 1 mv, this will not happen 
until 10‘^ Hz. 

As we show in a later section, there is no point in attempt- 
ing to observe hm waves below about 0.01 Hz. This could be done 
according to our previous argument if the instrument noise is about 

3 mv in the band 0.01 to 1 Hz. Our allowable resolution would be 

1 3 

roughly one part in 2 . We feel this is attainable with good 

electronic design. Consequently, it appears feasible to measure 
hm waves in Shuttle orbit with a search coil magnetometer. 
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B. Fluxqate Magnetometer 

If we must measure the absolute field with the same instru- 
ment as the hm waves, we have a much more difficult problem. In 
particular, we must use some device sensitive to the total field 
component along a given axis. Such devices as fluxgates or alkali 
vapor with biasing coils can accomplish this. Usually such in- 
struments operate as null detectors, feeding back a current through 
a biasing coil. The fundamental limitation in such instruments is 
the fluctuations of the feedback current that are observed by the 
null detector as field variations. An analysis of instrument noise, 
and hence feasibility for use in Shuttle orbit, depends on the 
particular instrument and strategy chosen. As an example, the 
UCLA fluxgate magnetometer developed for ground use resolves 65,536 
Y to one part in 2 , i.e., to 1/16 y (62 my). Instrument noise in 
the ULF band (10'^ to 1 Hz) is less than this resolution. We 
note, however, that the absolute accuracy is not better than 32 y. 

On a long time scale the instrument can drift by amounts comparable 
to this. Since the time scale of such drifts is days, it is of no 
consequence in ULF wave measurements. 

To demonstrate the feasibility of using such an instrument 
in Shuttle orbit, we briefly describe the design of the UCLA flux- 
gate magnetometer. As we will show, there is no reason why this 
design will not work with the same sensitivity in Shuttle orbit. 

A block diagram of a single axis is shown below. The mag- 
netometer includes a ring core sensor, second-harmonic fluxgate 
detector, analog feedback winding for field nulling, a level detector 
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an up-down counter, a di gi tal -to-anal og converter, and an offset 
winding for additional field nulling. The fluxgate detector acts 
as 
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a highly sensitive null detector of limited dynamic range. Pro- 
vided the field along the sensor does not exceed this dynamic 
range, the fluxgate detector generates a feedback current which 
completely nulls out the field along the sensor axis. If the 
field exceeds the dynamic range, the fluxgate detector output 
trips a level detector. This, in turn, enables an up-down counter 
which adds or subtracts successive counts from its current con- 
tents. The contents of this counter are continuously converted 
to a current proportional to the count. This current passes through 
an offset coil partially nulling the; field along the sensor. Through 
proper scaling each count corresponds exactly to half the dynamic 
range of the fluxgate detector. Thus, each time the detector 
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exceeds its dynamic range it is brought back to the center of this 
range by a change of one count in the offset field generator. 

In the UCLA ground magnetometer the offset field generator 
is scaled so that 11 bits encompass the complete range of pos- 
sible earth fields, i.e. 

2^^ bits = 65,536 y 

or 

2048 bits = 131,072 y 

Thus, the least significant bit (LSB) corresponds to 
LSB = 64 y 

The fluxgate detector is, in turn, scaled so that it has a dynamic 
range of ±64 y which is also encompassed by 11 bits. Thus 
2^^ bits = 64 y 

or 

LSB = 1/16 y 

The digital output of the instrument is obtained by appending 11 
zeroes to the binary count in the up-down counter and algebraical- 
ly adding the 11-bit output from the anal og-to-di gi tal converter. 

Instrument noise in this system is very difficult to evaluate 
theoretically. Past experience has shown that it depends heavily 
on the geometry of the circuits, the degree of isolation of the 
magnetometers' power supplies from other instruments (particularly 
digital), the degree of integration (used in circuit construction 
and improvements in solid state technology (particularly ADC's 
and DAC's). 
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To demonstrate the feasibility of this instrument we adopt 

the measured noise spectrum at the output of the fluxgate detector 

in the UCLA magnetometer as typical. This is found to be a 1/f 

-4 2 

spectrum with spectral density of 4 x lO y /Hz at 1 Hz. Accord- 
ing to our discussion of the natural spectrum given previously, 

an acceptable lower limit to the natural spectrum is a 1/f spec- 

-5 2 

trum of spectral density 2.8 x lo y /Hz at 1 Hz. These spectra 
are plotted below as a function of frequency. 
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Using the above estimates it is apparent that this fluxgate is 
usually adequate for frequencies below 0.1 Hz. However, above 
this frequency instrument noise exceeds the lower limit of the 
natural spectrum. However, in our earlier discussion of the 
search coil magnetometer we concluded the search coil would 
probably be inadequate below about 0.03 Hz. Thus, both search 
coil and fluxgate would be required to observe the entire ULF 
band. In a later section we will demonstrate that the lower 
the frequency the more likely it will be that spatial variations 
will be confused with temporal variations in Shuttle orbit. In 
particular, 0.01 Hz appears to be a lower limit of observable 
variations. Consequently, it appears that a fluxgate is not an 
appropriate sensor for hm wave measurements in Shuttle orbit. 

The origin of the noise in the fluxgate is not clear. If we 
assume it is due to electronics associated with the fluxgate detec- 
tor, we might be able to improve the overall noise by rescaling 

the instrument. For example, since this detector is scaled so 

-2 

that 64 Y - 10.24 v, we have a sensitivity of 16 x 10 v/y. Using 

2 2 

this to convert our observed noise spectrum in y /Wi to (v) /Hz, 
we obtain a noise spectral density at 1 Hz of 7 x 10“ (v) /Hz. 

This gives a two decade rms noise power of 1.8 mv. This is some- 
what better than we assumed possible for the search coil in a space- 
craft environment. 

If we reduce the dynamic range of the fluxgate detector by 
2 and assume that the electronic noise remains unchanged, we reduce 
the noise spectral density in y units by a factor of four. We 
note, however, that one additional bit must be added to the offset 
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field generator. Noise power from it will now be four times as 
important. In the present magnetometer design it appears that 
noise from the sensor, offset field generator and fluxgate detec- 
tor are comparable. Consequently, compromises such as those sug- 
gested above cannot significantly improve the performance. 

Thus, we conclude that simple modifications of the current 
fluxgate design are not likely to significantly change the in- 
strument noise. Consequently, the search coil remains the best 
choice for an hm wave sensor in Shuttle orbit. 
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Difficulties Associated with hm Wave Measurements at Shuttle Altitude 
In a preceding section we discussed various instruments with 
sufficient sensitivity to measure hm waves in space. If one of 
these instruments is placed on a spacecraft at Shuttle altitude, 
there will be a number of factors which limit its ability to detect 
such waves. In this section we consider only those factors which 
are important even if we had an ideal measurement platform. These 
are nearly all consequences of motion of the sensor around the 
earth. They include separation of the spatial and temporal features 
of the hm waves, effects of spatial gradients in the main field, 
crustal anomalies, ionospheric currents (Sq, equatorial electro- 
jet, auroral electrojet), and field aligned currents. 
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A. Separation of Spatial and Temporal Aspects of hm Waves 
Let us suppose that we are able to make perfect hm wave 
measurements on a platform in low altitude polar orbit. With a 
single moving platform it is impossible to distinguish between 
temporal changes of the wave and spatial changes. For example, 
suppose that a 100 sec period wave is localized in latitude to 
about 1000 km. At the velocity of a spacecraft in low altitude 
polar orbit (-vlO km/sec) it will take 100 sec to travel across 
the region of localization. In this time the wave amplitude at 
a fixed spatial point will have undergone one cycle. For the 
moving magnetometer it will be impossible to decide whether the 
observed change in ^ was due to temporal changes or spatial varia 
tions. 

To demonstrate the above, suppose the spatial variation in 
a given field component is Fourier analyzed. Then if the field 
is localized in latitude with scale x, the most important wave 
numberswill be those for which kx^ 2it, or X 'v x^, i.e., wave- 
lengths comparable to the scale. If the field component is quasi 
sinusoidal in time, we have 

(x, t) a exp [iC2iT/x)xJ exp [iwt] 
or 

(x, t) oc e [i2Tr(x/X + t/T)J , 
where w = 27r/T. But for a moving spacecraft x = vt, so 
B^(t) cs exp CIBttCv/X + 1/Tj t 

Clearly spatial variations in the field are seen by the moving 
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magnetometer as temporal changes. In particular the problem is 
most serious for wave periods T 'v x/v. 

For latitudinal localization of order 1000 km and satellite 
velocities 'v^lO km/sec, we have T 'v 100 sec. This is within the 
period of interest to us in studies of hm waves. 

In a recent study of a stormtime Pc 5 wave event, Lanzerotti 
^ £l_. (1 974) examined the latitude dependence of a hm wave of 
400 sec period. Their analysis suggests that the wave had a 
spatial scale of 'vlOOO km (10° of latitude) centered about 62° 
magnetic latitude. From our preceding analysis it is clear that 
time and space variations would have been badly mixed in a Shuttle 
observation of this wave event. 
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B. Effects of Motion Through the Earth’s Dipole Field 
As a polar orbiting satellite passes around the earth, the 
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direction and magnitude of the field at the satellite are con- 
tinually changing. Because the earth's field is so large, the 
rate of change of field due to satellite motion dominates any 
changes due to hm waves. To demonstrate this, consider the 
fol lowing. 

In spherical coordinates the dipole field is given by 
B^ = -( 2 Bjj/L^) cos 0 
Be = sin e 



The output of a search coil sensor aligned in the radial di- 
rection will be proportional to dB^/dt. 

C(d/dt)B^] = [(a/ 3 t)B^ + . ^)B^] 

We assume there is no time variation at a fixed point, so (3/3t)B^ 
= 0. Expanding • ^)B^ in spherical coordinates gives 

(d/dt)B^ = [(;^ • = v^/(l )(a/ 3 r)B^ +(^ 9 /^) (3/30)B^ 

v^/r sin e ( 3 / 34 >)B^ + [B^/CDCr)] [v^ (3/30)(l) - Vq 

( 3 / 3 r)(r)] + [B^/(l)(r sin- e)] [v^(3(l )/3(|)) ^ v^ 

(3(r sin e)/ 3 r)] 

But B. = 0 and we assume for simplicity a perfect circular, polar 
9 

orbit with v^ = v. = 0. Thus 
r (|» 
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(d/dt)B^ = (Vq/t) C a/3e)B^ + Bg/r [-v^J = (v^/r)(2B^ sin 8/L^) 
-(VqBq/t) (d/dt)B^ =(Vg/r)(-2Bg) - (v^B^/r) = -3VgBg/r 

Finally, 

(d/dt)B^ = (3Vg/r) (B^/L^) sin e = (3Vg/Re ) (B^/ (L )^ ) sin e 

For a near-earth circular orbit L 'u 1 and e = fit, where 

= v./Re. Thus 

0 

[(d/dt)B^]L^^ = C3^S®o^ 

The output of the sensor is a sine wave with amplitude 

3 fl 5 Bjj = 50 TT y/sec . - =■ ' 

This is approximately 160 y/sec. 

For comparison, consider a hm wave of 106 sec period and 
2.5 Y amplitude. B = B^ sin wt so, 3B/9t cos wt. 

The output of the sensor is again a sine wave with amplitude 

o)B„. 

o)By = (27T/T)By = 2tt/100 • (2.5) 'v 0.157 y/sec 
The ratio of these signals is 

l/(signal wave/signal dipole) = 1 (2TrBy/T)/50Tr = 25 T/Bj^ 1000. 
Clearly even a large amplitude hm wave has an effect a 
thousand times smaller than the effects of motion through the 
earth's field. 
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Attempts to remove the earth's field using field models and 
information on the sattelite orbit can introduce serious problems 
in the data from a satellite in shuttle orbit. An example shown 
in Figure 1 of one such attempt made by Cain et al , 1967, 
illustrates the problem. Residuals left after subtracting the total 
field of a model from the observations show large amplitude oscil- 
lations in the field magnitude. The correlation of these with 
satellite perigee makes the results very suspicious. Possible 
sources of error include higher order terms in the field model 
and small errors in the determination of the satellite orbit. 

For a vector field instrument the problem would be more serious 
since errors in atttitude would enter as well. Clearly it would 
not be easy to determine whether these oscillations were waves 
of natural orgin or not. 
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C. Effects of Crustal Anomalies 

The earth's field has higher order terms than a dipole, 
which will also affect the measurements made by an hm wave sensor 
at Shuttle altitude. Furthermore, close to the earth's surface 
anomalies caused by variations in the composition of the crust 
become important. Since the magnetic field due to these anomalies 
satisfies Laplace's equation, it is possible to calculate the 
field at Shuttle altitude if it is known at the earth's surface. 
Unfortunately, very little is currently known about the spatial 
distribution of anomalies of scale important at Shuttle altitude. 

The reason for this lack of knowledge is that the process 
of upward continuation of the surface field acts like a low pass 
filter. Roughly, one can say that only features of scale comparable 
to the distance from the source to the observation point can be 
observed. At Shuttle altitude this is 500-1000 km. This distance 
is so large that, to date, there are no surface magnetic maps 
showing features of this scale. The reason is simply that it 
takes so long to survey a region of this size that temporal varia- 
tions in the surface field dominate the observations. With the 
current sparse distribution of fixed observatories it is impossible 
to remove these temporal effects from the survey observations. 

There has been one attempt to generate a map of crustal 
anomalies using satellite data, however. Zietz e^ aj,. (1970) used 
data from the USSR satellite Cosmos 49 to contour anomalies ob- 
served between 261 and 488 km altitude over the United States. 

Their results, shown in Figure 2 suggest that a significant problem 
exists for an hm wave sensor at Shuttle altitude. It is evident 
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from this figure that a satellite pass over the U.S. would observe 
quasi peri odi c variations of order 50 y in distances of order 1000 
km. At 10 km/sec satellite velocity these would appear to be hm 
waves of 50 y amplitude and 100 sec period. 
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; D. Effects of Ionospheric and Field Aligned Currents 

For a polar orbiting satellite at a few hundred km altitude, 
ionospheric currents will cause significant variations in the 
field. Distinct currents which must be considered include the 
solar quiet day variation, Sq; the polar cap currents due to mag- 
netospheric convection, Sqp, the equatorial electrojet and the 
auroral electrojets. All of these current systems have con- 
siderable spatial structure and are time varying as well. This is 
particularly true of the auroral electrojet with a latitudinal ex- 
tent of a few hundred km and time variations on the minute time 
scale . 

Effects of the Sq current system at the earth's surface 100 
km below the source are of order 50 y occasions this current 
system has spatial structure, with scale tens of degrees (several 
thousand kilometers) (Schieldge, 1974). Effects of the equatorial 
electrojet are even larger at the earth's surface and have been 
studied by satellite magnetometers at Shuttle altitude (Cain and 
Sweeney, 1973). Figures 3 and 4 taken from this paper illustrate 
the nature of the observed variations in total field. Changes as 
large as 20 y seen in less than 5° of latitude as the satel- 
lite passes over the magnetic equator. This corresponds to 500 
km or 50 sec, i.e., a rate of change of order 0.4 y/sec. Clearly, 
the variations appear to be quasiperiodic and would be indistinguish- 
able from hm waves of long period. Comparison of different passes 
over the equator show that this phenomenon is highly variable and 
probably impossible to predict and remove from the data. 

f«rturbations due to the auroral electrojet are much larger 
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than those of the equatorial electrojet at Shuttle altitude. 
Figure 5 taken from a paper by Lange! and Cain (1968) shows the 
effects in the total field. For the pass shown there was a 
change of 500 y within 10® of latitude. 

Passage of the satellite through sheets of field-aligned 
current flowing into the auroral electrojet will also cause per- 
turbations of large amplitude. In a recent paper by Armstrong 
e^ aj_. (1974) vector field observations made by a fluxgate mag- 
netometer on the TRIAD satellite as it flew through an auroral 
arc were reported. Figure 6 taken from this paper shows that 
in approximately 20 sec the east component of the field decreased 
by almost 400 y. Superimposed on this decrease were small-scale 
perturbations of many tens of y amplitude, presumably due to 
fine scale structure in the field-aligned currents. 
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E . Summary 

In the foregoing discussion we have outlined a number of 
sources of magnetic field variation which would be recorded by a 
satellite in a low altitude polar orbit. All of these sources 
would cause perturbations that would have to be removed from the 
data in order to make an unambiguous identification of an hm wave 
in the data. As shown above, most of these sources are unpredict- 
able, have relatively small spatial scale, significant time varia- 
tions and magnitude much larger than the waves of interest. As- 
suming that a strategy is devised for removing these effects, one 
is still faced with the problem of separating the spatial aspects 
of the perturbations from the temporal. 


32 


Difficulties Associated with Measurements on Space Shuttle 
A . AC Interference Created b.y Spacecraft 
The difficulty most likely to be encountered on the Space 
Shuttle is low-frequency ac noise generated by the spacecraft it- 
self. The Shuttle will be an extremely complex device constructed 
with magnetic materials, many of which will be moved as sensors 
are pointed or booms are deployed and retracted. There will be 
a large variety of electric motors used to actuate devices as well 
as numerous complex electrical circuits. It seems unlikely that 
much attention can be paid during spacecraft construction to mini- 
mizing the effects of these different sources of ac interference. 

An example of such circumstances of high levels of magnetic 
noise is that of the UCLA fluxgate magnetometer on ATS 1. This 
spacecraft was not designed for a low noise environment and the 
magnetometer was attached to the spacecraft itself. Since the 
spacecraft electronics were located only a few feet from the mag- 
netometer, changes in the state of various subassemblies were de- 
tected as changes in ambient field. Dc offsets of 50 y were 
regularly recorded in a sensor parallel to the satellite spin axis 
as one or another system was turned on and off. More serious from 
the point of view of wave measurements was an 82-second period, 
sawtooth variation in the dc field. Although the amplitude of 
thi s 1 variati on was only about 3 y pp, spectral analysis of waveform 

i 

produced a harmonic spectrum spanninig most of the band of hydro- 
magnetic waves. Finite resolution of the spectral analysis, 
smootiting of the spectral estimates, etc., spread the power in 
these harmonics producing a background spectrum above which a natural 
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signal must rise in order to be observed. 

A problem similar to the above was observed at higher fre- 
quencies in the transverse field sensor on ATS 1. In this case 
a ramp of period equal to the basic telemetry sequence (5.12 sec) 
was present with amplitude of order 5 y- Much of the time this 
signal dominated the waveform of the magnetometer making it im- 
possible to identify visually the occurrence of hm waves in the 
Pc 1 band. 

Examples of these problems are shown in Figures 7 to 9. In 
Figure 7 a sawtooth of several gamma amplitude in the Y component 
is phased with the vertical lines corresponding to the start of 
telemetry sequences. In the spectra of Figure 8 at least five 
harmonics of this 5.12-sec periodic waveform are evident. Figure 7 
also shows the longer period ramp in the Z component (see step in 
Z at second vertical line). At least 18 harmonics are present 
in the spectrum, shown in Figure 9. Finally, Figure 7 shows the 
transient effects of a dc offset in the field at the eleventh 
vertical line from the right edge. 

These examples clearly demonstrate that a magnetometer lo- 
cated close to an improperly designed spacecraft will be influenced 
by a variety of noise sources. We believe it is extremely unlikely 
that Shuttle can be constructed to avoid such problems. Conse- 
quently, we conclude that an hm wave sensor cannot be mounted on 
the Shuttle itself. In the following section we consider the 
possibility of using a long boom. 
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B. Problems Associated with Long Booms 

From the foregoing discussion, it is clear that the hm wave 
sensor should be located at some distance from the Shuttle itself. 
Long booms are one means of accomplishing this, provided it is 
physically possible to mount the sensor on the boom. Search coil 
sensitivity is roughly proportional to weight. Typical coils used 
for ground measurements weigh 100 lbs. apiece. A 300 pound sensor 
array mounted near the end of a 50 meter boom would generate a 
considerable moment of inertia when any attempt is made to change 
the spacecraft attitude. 

If we assume that sensor weight is not a consideration, we 
must sill consider the importance of boom vibrations. For example, 
let us examine the situation of a sensor oriented at an angle 
0 = 0 Q + A 6 to the earth's field, B^. The component of field along 
the sensor is then 

Bj = Bq cos (0q + A0) 

Expanding the cosine function gives 

Be = (cos 0„ cos A0 - sin 0 ^ sin A0) 

5 0 0 0 

Since A 0 is a small angle, we have 

Be 3 cos 0^ (1) - B- sin 0^ • A0 
S o o'' 0 0 

The magnitude of the perturbation due to boom oscillations A0 is 
maximum for a sensor perpendicular to the field, 0 ^ = 90®. 

AB 3 (90®) = C2it/ 360 • 50,000 y] • A 0 (deg) 
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or 


ABg (90°) = 870 Y * A0 (deg) 

This is roughly 14 y per minute of arc. 

The foregoing is a serious problem which has limited the 
ability of previous magnetometers to sense hm waves close to the 
earth. For example, the UCLA fluxgate magnetometer on OGO 5 was 
mounted on a 20-foot boom. As the satellite approached perigee, 
frequent changes in attitude were made to keep the antennas pointed 
earthward and the solar cell arrays perpendicular to the sun. 
Impulses caused by firing of attitude control jets caused damped 
sinusoids of 0.3 and 3 Hz to be generated. As many as 20 cycles 
of 5-y amplitude oscillations were recorded in a 1000-y field. 

These imply boom oscillations of order 0.25°. From our calculations 
above, this would give boom oscillations of 200 y at Shuttle. al- 
titude. Clearly, the booms used on Shuttle would have to be ex- 
tremely rigid if attitude control maneuvers are continuously car- 
ried out. 

It should be pointed out that boom vibrations generate high- 
ly coherent, monochromatic signals. In addition, since the vi- 
bration cannot alter the ambient field, the field magnitude must 
remain constant. Thus, the oscillations will have no compres- 
sional component. As a consequence of these facts, their mag- 
netic signature is easily recogni zed? wi th spectral analysis. Their 
importance arises from the possibility of saturating search coil 
amplitudes or obscuring weaker signals in waveform plots. 
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C. Problems Associated with the Determination and Control 
of Spacecraft Attitude 

Difficulties may arise in hm wave measurements because of 
uncertainties in Shuttle attitude. Present plans are to deter- 
mine this to about 0.5® accuracy. Using the calculations of the 
preceding section, this indicates a possible error in some com- 
ponent of the field as large as 400 y. In addition, we note that 
the booms to be used on Shuttle will be more than 100 feet in 
length. Maintaining sensor alignment at the end of such a boom 
seems likely to be quite difficult. It may be necessary to have 
a separate attitude determination system at the end of the boom 
to determine sensor alignment with sufficient accuracy. 

Errors in sensor attitude are most serious for a device 
measuring the vector field rather than its derivative. The problem 
arises when we try to use two or more satellites in the same or- 
bit as a gradiometer. Subtraction of the vector field observa- 
tions made at the two different satellites at a given time can 
produce a large difference vector. This difference will include 
the random errors due to attitude determination which would add to 
a difference vector as large as 1000 y. 

Changes in spacecraft attitude may be serious for a search 
coil if attitude control maneuvers are carried out in discrete 
steps rather than continuously. For example, suppose the rate of 
rotaJtion is 1° per second (6 min rotation period). From above, 
the rate of change of field could be as high as 870 y/sec, which 
wouldiisaturate any search coil system. 
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D. Summary 

Measurement of hm waves on the space Shuttle itself will be 
quite difficult for several reasons. First, the Shuttle is likely 
to be a source of low-frequency noise that will contaminate the 
observations. This noise could originate from a wide variety of 
electrical circuits, motors, and slowly moving magnetic materials. 
In order to reduce the effects of this noise, the hm wave sensor 
could be placed on a long boom. However, if the sensor consists 
of a highly sensitive search coil, it will be very heavy. This 
could adversely affect the moment of inertia of the spacecraft. 

Even if this is not a problem, vibrations of the boom as the 
attitude of the spacecraft changes could produce extremely large 
field oscillations that would either saturate the measuring 
system or swamp any natural signal in the data. Assuming that pro- 
per boom construction can damp the oscillations, attitude control 
maneuvers can produce rates of change of field that would saturate 
the system. Errors in attitude are another problem, particularly 
when more than one satellite in the same orbit is being compared. 
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Difficulties Associated with Measurements 
on Tethered and Free Subsatellites 

In the preceding section we showed that it is difficult to 
measure hm waves on the Shuttle itself. Two ways in which this 
problem might be solved are to use either a tethered subsatellite 
or a free subsatellite. As we will show, there are difficulties 
associated with both of these. 

A. Tethered Subsatellite 

A tethered subsatellite appears to have a number of advan- 
tages for hm wave measurements. First, the subsatellite can be 
located far enough from the Shuttle to minimize low frequency 
noise. Second, power to the subsatellite and data communications 
can be through the tether to the Shuttle. Third, the attitude 
of the subsatellite can be maintained fixed in inertial space. 

In order to realize the foregoing advantages, there are 
several problems which must be considered. First, the tether 
would have to be attached to the Shuttle along an axis which re- 
mains fixed in inertail space. If this is not the case, any at- 
titude control maneuver would set the combined system of Shuttle 
and subsatellite into rotation. This could be done, for example, 
by orienting the longitudinal axis of the Shuttle transverse to 
the Shuttle velocity and tethering the subsatellite to the nose 
of the Shuttle. As the Shuttle traveled around the earth, it could 
roll so that the top (or bottom) of the Shuttle was always along 
the local vertical. 

The foregoing configuration could cause another problem to 
arise', however. Because of the size and shape of the Shuttle, 
drag from the atmosphere will be greater on it than on the sub- 



satellite. The unbalanced forces on the two ends of the tethered 
system constitute a torque that would set the combined system 
into rotation about the center of mass. How significant this would 
be depends on the altitude of the orbit and the shapes of the 
tethered Shuttle and subsatellite. 

Another problem which must be considered is how to keep the 
subsatellite fixed in inertial space. Clearly, spin stabilization 
is not possible because of the large field changes this would 
cause in the hm wave sensor. If the subsatellite carries a heavy 
array of three orthogonal search coils the subsatellite might 
have a tendency to tumble as a consequence of forces due to its 
own attitude control system or impulses transmitted to it via 
the tether. It seems likely that gyros may be required to sta- 
bilize the subsatellite. Finally, we note that there must be a very 
accurate system for determining the attitude on this subsatellite 
to allow reduction of observations from multiple satellites to a 
single coordinate system. 
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B. Free Subsatellites 

Free subsatellites are another alternative which might 
make it possible to observe hm waves in Shuttle orbit. Clearly, 
the problems of determining and maintaining fixed inertial attitude 
of a free subsatellite is the same as for a tethered subsatellite. 
In addition, however, there are new problems of power and com- 
munications. For a short-lived mission the power problem does not 
appear to be too serious, but it should be kept in mind that the 
subsatellite will be in darkness half the time. Communications 
are more difficult since it cannot be assumed that ground stations 
will always be under the subsatellite. Possible solutions in- 
clude on-board data storage, communications via synchronous sat- 
ellite and communications directly to or through other subsatel- 
lites with the Shuttle. 

Another problem with a free subsatellite is maintenance of 
its position in orbit relative to the Shuttle and other free 
subsatellites. For short missions this would probably not be 
serious, but if a number of subsatellites are left in orbit after 
Shuttle returns to earth, the problem would become important. 
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C . Summary 

Both tethered and free subsatellites provide a better op- 
portunity for hm wave measurements than Shuttle itself. Both 
types of subsatellites require a method of precisely determining 
and maintaining the attitude of the subsatellite in inertial space. 
The interaction of the Shuttle with the subsatellite via the te- 
ther is a problem of some concern. The tether does, however, pro- 
vide a simple means of providing power and communications to the 
subsatellite. A free subsatellite requires a source of power, 
means for controlling its orbital position, and communications 
with the ground, considerably increasing the cost and complexity 
of the free subsatellite as compared with a tethered subsatellite. 
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Magnetometer Array for Separating Spatial 

and 

Temporal Effects in Hm Wave Measurements at Shuttle Orbit 
We have previously dsicussed a number of sources of magnetic 
field which would be observed by an hm wave sensor in Shuttle or- 
bit. It is apparent that satellite motion through any spatial 
changes in the magnetic field would be recorded as a time variation. 
Furthermore, the absence of time variations at a particular ground 
station beneath the satellite could not be taken as evidence that 
the satellite signal is due to purely spatial effects because of 
the possibility that wave phenomena are spatially localized. Be- 
cause of these facts, we must devise a strategy for separating 
spatial and temporal effects. As we show below, this requires an 
array of magnetometers if the results are to be unambiguous. 

A. Observations with a Single Observatory 
To begin, let us consider the simpler problem of separating 
temporal and spatial effects in ground observations. If we have 
a single observatory, we can only measure the vector field as a 
function of time. If observations are made over a sufficiently 
long time, the mean value of the observations will be due to spa- 
tial effects, i.e., local sources of magnetic field. The fluctua- 
tions about this mean will have various causes including propagat- 
ing waves, standing waves, and time varying noise. We may Fourier 
transform the observations obtaining the frequency spectrum of 
the fluctuations. If the Fourier spectrum were entirely due to 
waves and we knew their dispersion relation, we could calculate 
the spatial distribution of the waves as well. This follows from 
the fact that the waves satisfy the wave equation which enables 
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one to write for any field component 

In the above k(u)) is the dispersion relation for the particular 
wave phenomenon and A(w) is the Fourier transform of the observa- 
tions at a particular location. 

Since our major experimental goal is to define the dispersion 
relation, we must invert the above procedure. In other words, we 
wish to use observations of f(x, t) to define k(w). Thus, we must 
make observations at more than one point as a function of time. 
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B. Observations with Two Observatories 

Suppose the signal consists of a single frequency component, 
o)q. Then 


Note that this implies that only a single wavelength is involved. 
From the observations at a location X, we can determine the fre- 
quency Wq. Then from observations at a second location we can 
determine the difference in phase, i.e., 
f (x^ , t) = cos (kx-j - wt) 

f(x2» t) = cos (kx2 - o)t) 

A(J) = k(x2 - ) = kAx 

Then k = A(j>/Ax and = u/k. For this simple case, two ob- 

servatories are sufficient to determine the wave number and phase 
velocity. If nature were sufficiently kind as to provide repeated 
monochromatic examples of this wave phenomenon over a wide range of 
frequencies, we could map out the dispersion relation. 

In many natural situations we may have A<() = 0 , i.e., standing 
waves. If this happens, we learn nothing about the dispersion re- 
lation. More frequently the signal consists of a broad spectrum of 
waves and no constant phase difference exists. Again, we are un- 
able to determine the dispersion relation. 
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C . Observations with a "Dense" Array of Observatories 
In the preceding discussion, our ability to determine the 
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dispersion relation depended' on an assumption about the shape of 
the signal waveform. In order to solve the general problem we 
must make independent measurements of both the time variations 
and the spatial variations of the signal. To do this, we must 
have an array of observatories. 

Let us assume for the present that we can continuously mea- 
sure both the time and space variations of our signal with infinite 
resolution. In this case, we can evaluate the double Fourier 
transform for the observations 



This transform may be displayed as contour maps of the real and 



In such a display any fluctuations which satisfy the wave equation 
will have f(k', w) organized along trajectories corresponding to 
the dispersion relation. This follows from the fact that for a 
given only those values of k which satisfy the dispersion 
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relation are allowed. Note it is this fact which enabled us to 
express f(x, t) as a single integral when we discussed the spatial 
distribution earlier. 

The plot of f(k, w) in the w, k plane is an extremely useful 
tool because it allows us to distinguish between propagating waves, 
standing waves, noise and purely spatial variations. First, since 
purely spatial variations are time independent by definition, 
they appear at zero frequency. Because they are unrelated to the 
wave phenomena present in the signal, they will appear as a singu- 
larity along the k axis. Second, noise will not satisfy the wave 
equation and should not be organized along trajectories. Third, 
monochromatic waves appear as a pair of impulses in opposite qua- 
drants. The phase velocity of this wave is given by the slope of 
a line from the origin to the locations of these two impulses. 
Finally, since a standing wave is the superposition of two waves 
propagating in opposite directions, it will be represented by 
two pairs of impulses symmetrically located in all four quadrants. 
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D . Limitations Resulting from Finite Extent and Separation 
of Observatories in Array 

In the discussion above, we assumed that we could measure 
the spatial distribution of the waves with infinite resolution 
over all space and all time. In a typical situation both the 
length of the array and the time duration of the records are of 
limited extent. Furthermore, the data are usually sampled at 
discrete intervals. These effects of truncation and discrete 
sampling can be quite significant, particularly in the case of 
the spatial variable where each additional observatory represents 
a considerable increase in cost. 

The effects of truncating a time or space series in the 
transform domain are easily calculated. Truncation is accomplished 
by multiplying the series by a "rectangular" function. Using 
the convolution theorem, multiplication in one domain is equiva- 
lent to convolution in the other. Thus, we convolve the sin X/X 
transform of the rectangle function with the transform of the 
original data. For a rectangle of length L the transform will have 
a width of order 1/L. As a result of the convol ution ,detai 1 - in 
the spectrum of width less than 1/L is lost. In particular, we 
cannot determine any information about frequencies less than 1/L, 
i.e., wavelengths longer than the length of the array. 

Effects of sampling are calculated in a similar way. Sampling 
is accomplished by multiplying the series by a Dirac comb of 
delta functions spaced AL apart. The transform of the Dirac comb 
is also a Dirac comb with reciprocal spacing, 1/AL. Convolving 
this transform with that of the original series replicates the 
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spectrum an infinite number of times. Problems arise if the spec- 
trum of the original data extends more than half way to the first 
delta function of the reciprocal Dirac.comb, i.e., beyond 1/2 aL. 

When this is the case, convolution folds the spectrum on top of 
itself and it is impossible to determine whether the original 
signal is above or below the folding frequency. In the original 
domain this process of folding corresponds to sampling the ori- 
ginal waveform at intervals greater than half its shortest wave- 
length. 

If observations are made with a fixed ground array, the main 
limitation is the number of observatories which is available. 

Usually we can sample with as much time resolution and for as 
long a time as necessary. Consequently, the frequency resolu- 
tion and folding frequency that result from truncation and sampling 
in the time domain are not important. In the wave number domain, 
however, serious problems arise. For example, suppose 10 observa- 
tories are used, then only five wa vel engthsccan be measured. In 
many cases this would be adequate to determine considerable in- 
formation about the dispersion relation. Considerable improvement 
could be achieved if one adopts the strategy of removing every 
other station and using them to double the lengths of the array. 

Over a long interval of time this procedure would define the 
phase velocities of a large number of wavelengths. 

The foregoing procedure provides a possible means for sep- 
arating temporal and spatial effects in Shuttle orbit. Suppose 
that the Shuttle places a number of subsatellites into the same 
orbit each separated by a variable distance. Data from this Shuttle 
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array could be double Fourier transformed and reduced to a sta- 
tionary coordinate system using the shift theorem of Fourier 
analysis. This theorem states that the Fourier transform of a 
function shifted by an amount a, i.e., f(x - a) is given by 

where 

For the Shuttle array, a = V^t, where = satellite velocity. 

The need to define the time variation by a moving array im- 
poses a new restriction not discussed above. Since the array is 
of finite extent the time variation in a given spatial interval 
can only be measured while the array is passing by. This fact is 
illustrated below where we show the location of the array 

as a function of time by diagonal lines in the x-t plane 



From the diagram it is clear that no time variations can be mea- 
sure*l‘:iat a given spatial location af,ter the last observatory has 
passed overhead. Note that the longer time, T, one attempts to 
measure, the shorter is the spatial interval 0 over which this can 
be done. If we optimize both D and T by maximizing the area of 
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the rectangle, then D =L/2 and T = L/2 V^. As an example, suppose 

4 

the length of the array is 10 km and 10 observatories are used. 

Then from above the effective length is 5000 km and T = 500 seconds. 
Note that at any one time only half the observatories are con- 
tributing to the definition of the function of x and t. Because of 
this, only half as many wavelengths may be defined, as was the case 
for the fixed array. Thus, 10 moving observatories will define 
only two wavelengths. 

Provided it is possible to progressively increase the sepa- 
ration of the observatories, it would still be possible to define 
the dispersion relation of a frequently occurring wave phenomenon. 
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Summary of Feasibility Study 

The fundamental question we have examined in this report is 
whether it is feasible to measure hm waves on space Shuttle. The 
answer appears to be that it is just possible technically to make 
such a measurement. However, for a number of reasons, the results 
would be very difficult to interpret. 

The basic technical problem is that hm waves are extremely 
weak in comparison with the earth's field at Shuttle orbit. Any 
instrument designed to measure these waves must simultaneously 
have a wide dynamic range and high resolution. For example, a 
fluxgate magnetometer would have to have a resolution of order 
0.01 Y dynamic range of +50,000 y, or 1 part in 10^. 

For a search coil on a non-spinning spacecraft, the require- 
ments are less severe, but still significant. The primary re- 
quirement is that the induced emf due to satellite motion through 
the earth's field does not saturate the system amplifier. This 
requires that we be able to measure a signal of amplitude 
due to a wave superimposed on a signal as large as 50 iry/sec, the 

maximum dB/dt due to satellite motion in the dipole field. Thus, 

-4 

we must resolve to better than 2irB^/50ir T^ = 4 x 10 for a 0.01 y 
wave at 1 Hz. 

If the spacecraft is subject to frequent changes in attitude, 
the measurement becomes more difficult for either type of instru- 
ment. For a search coil we must still be able to resolve B,. 
compared with fly.ot®o* better than(2Tr6 /T^ J/2 ttBq/T^q^ = 

(Vot/Twave> 

If we assume that a typical attitude control maneuver is 
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carried out at a rate corresponding to a full roll In 500 sec, 
then resolution = (500/1) (10"^/5 x 10^) = 1 x 10'^ for the 0.01 
Y wave at 1 Hz. 

While these requirements are severe, existing instruments 
can probably be modified to meet them. Assuming this is done, 
we still believe it is extremely unlikely that hm wave measure- 
ments can be made on Shuttle itself. this is because the Shuttle 
observatory will be. so large and complex that there will be abun- 
dant sources of DC and AC field to contaminate the observations. 
Since the guiding philosophy of Shuttle is to be a low cost, general 
purpose facility there is little possibility that it will be con- 
structed to minimize these sources of noise. 

Long booms are one possible way to reduce the effects of 
Shuttle noise sources to a tolerable level. Such booms introduce 
a number of problems, however. The most serious problem is boom 
vibrations. A long boom is a resonant structure which will oscil- 
late in response to each change in attitude of the Shuttle vehicle. 
Effects of these vibrations could be much more serious than those 
due to both satellite motion through the dipole field and attitude 
rotation, depending on their amplitude and period. For example, 
if the boom vibrations were only 10 minutes of arc in 10 sec, the 
amplitude of the corresponding field oscillations would be about 
150 Y* This is enormous compared with any naturally occurring 
signal of this period. Spectral analysis would be required in 
order to extract any natural signals from the resulting waveform. 

A more attractive solution to the problems imposed by noise 
and attitude control on the Shuttle is to place the hm wave sensor 


53 


on a subsatel 1 1 te . If the subsatellite were tethered to the 
space Shuttle, power, ground communications, and orbital control 
could be supplied by the Shuttle, simplifying construction of the 
subsatellite. For a free subsatellite these would have to be built 
into the subsatellite in addition to other requirements. The most 
important of these requirements for magnetic field measurements 
would be inertial stabilization and a very accurate means of de- 
termining sensor attitude. If the instrument measures the ab- 
solute field rather than its derivative, very negligible errors 
in attitude (i.e., of order 1/10®) would cause large errors in 
the field components (as large as 87 y)* If measurements on 
two satellites were compared with attitude determined no more ac- 
curately than the 1/2° planned for the Shuttle itself, differences 
in measured field as large as 1000 y could result. Again spectral 
analysis would be required to eliminate these differences in the 
measured signals. We note that the requirement that the subsatel- 
lites be inertially stabilized may be hard to realize, particularly 
for a tethered subsatellite because of forces transmitted to it 
from the Shuttle. 

Let us now assume that a suitable designed magnetometer is 
placed in an inertially stabilized, noise-free subsatellite. 

There are numerous natural sources of magnetic field which will 
be observed by the instrument as time variations with amplitude 
greater than thse of hm waves. These include motion of the satel- 
lite over crustal anomalies in the earth's field, motion over the 
ionospheric Sq current system, motion over the equatorial and 
auroral electrojets, and passage through sheets of field aligned 
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currents. External to the satellite are the magnetopause, ring 
and tail currents. Effects due to the latter are not likely to 
be localized in latitude but they do change on a time scale 
comparable to the satellite's orbital period. The electrojets 
can be highly structured and passage over them would be recorded 
as large amplitude waves of short duration (several hundred gamma 
i n tend of seconds ) . 

None of the foregoing effects can easily be removed from 
the data of a single satellite. Most of them are due to phenomena 
which are highly variable on a day-to-day basis. Even the crustal 
anomalies, which are essentially constant in time, would require 
very long intervals of recording, perfectly circular orbits and 
synchronized orbital periods to remove from the data. If we limit 
our observations to the =:45® latitude between electrojets, we can 
only observe waves of 500 sec period or less since this is the 
time the satellite takes to traverse this region. For a typical 
Shuttle orbit of 500 km altitude, any observed time variation with 
period longer than 500 km/10 km/sec = 50 sec may be due to crustal 
anomalies. For periods shorter than this, it is more likely they 
are real time variations since far from a magnetic source only 
wavelengths comparable to or larger than the distance of observa- 
tion can be observed. 

Let us now assume that we have successfully measured some- 

i 

thing between 10 and 500 cycles of what are probably hm waves of 
period less than 50 sec^ These measurements by themselves are 
of little value, but if they are combined with ground data, then 
it becomes possible to make inferences about wave propagation. 
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source regions, etc. As an example, suppose nearly the same fre- 
quency is seen at both ground and satellite. A cycle-by-cycle 
comparison of wave polarization between satellite and ground start- 
ing when the satellite passed poleward over the station might re- 
veal a progressive change in phase due to wave propagation, or 
a systematic change in amplitude due to a standing wave. 

Since one observation is made above the ionosphere and one 
below, it is quite possible that some of the observed differences 
are a result of propagation through the local ionosphere. In 
order to remove possible ionospheric effects from the observations, 
one needs at least two satellites in the same orbit. Intercom- 
parisons between these and a ground station would reveal additional 
informati on . 

In general, most wave phenomena consist of a broad band of 
frequencies and wavelengths. In such cases a comparison of the 
waveforms at two locations will not be particularly meaningful. 

This is especially true if the medium is dispersive. In this 
case the only recourse is to establish an array of observatories. 

The largest wavelength which can be studied corresponds to the 
length of the array; the shortest wavelength to twice the array 
spacing. Similarly, the lowest and highest frequencies are limited 
by the duration and sampling of the time record at each station 
in the array. 

If we assume that data are available from such an array, we 
can best study wave phenomena by Fourier transforming the data to 
frequ'ency-wavenumber space, i.e., the w-k plane. A wave packet 
pi^opaga,t.ing through a dispersive medium will have a Fourier transform 
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organized along a curved trajectory in the w-k plane. A standing 
wave would consists of a symmetric pair of such trajectories cor- 
responding to pairs of waves at each frequency propagating in 
opposite directions. Truly spatial variations will appear as a 
singularity function along the k axis. 

If we wish to make an unambiguous determination of the phase 
velocity of hm waves, we must carry out such a procedure. If 
this were our only goal, it would be sufficient to use a fixed 
ground array. But if we wish to remove effects of satellite mo- 
tion over spatial variations from temporal effects, the array 
must be in Shuttle orbit. As we mentioned above, this is a matter 
of necessity for wave periods longer than about 50 sec. Unfor- 
tunately, a moving array of finite extent is handicapped by the fact 
that it is passing over a given spatial interval for a limited time. 
Since the double Fourier transform effectively requires equal 
length time records at each spatial point, we are only able to use 
data from half the spatial array. Thus, definition of only four 
wavelengths would require 4 x 4 = 16 observatories in the Shuttle 
array. 

It seems extremely unlikely that inertially stabilized sub- 
satellites capable of obtaining and maintaining their positions 
in a spaced array can be made and launched cheaply. Consequently, 
we feel a Shuttle array is unfeasible. The only alternative 
appears to be a detailed mapping of the vector magnetic field 
with a* single satellite in Shuttle orbit. Such a map would require 
an extensive interval of data acquisition and reduction to produce. 
Given such a map, knowledge of satellite position and attitude 



would allow one 
servations of a 


to remove purely spatial 
given orbit. 
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variations from the ob- 
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1. An hm wave sensor is technically feasible for Shuttle orbit. 

2. Existing sensors are inadequate in terms of resolution, 
dynamic range, and frequency response, but can probably be modi- 
fied to make the necessary measurements. 

3. It would be impossible to mount the sensor on the Shuttle 
itself because of high levels of magnetic noise. 

4. A long boom attached to the Shuttle is unlikely to be 
very helpful unless attitude control maneuvers are made continu- 
ously rather than in discrete steps. High rotation rates or 
boom ! vibrations are likely to mask any natural signals. 

5. A tethered subsatellite appears to be an inexpensive 
way of removing the hm sensor from the influence of the Shuttle 
provided it can be inertially stabilized and will not be influenced 
through the tether by attitude control maneuvers of the Shuttle. 

6. A free subsatellite that can be positioned as well as 
stabilized is a better location for an hm wave sensor. 

7. Hm wave measurements can probably only be made in the 
s:45® of latitude between the highly structured and' unpredi ctabl e 
equatorial and auroral electrojets. 

8. Because of magnetic effects due to crustal anomalies, 
wave periods longer than ==50 sec would be difficult to study. 

9. Studies of long period waves would require either an 
array of sensors in Shuttle orbit or a long-term mapping of the 
crustal anomalies. 

10. Effective wave studies would require at least two variably 
spaced sensors in Shuttle orbit and one ground station. 

11. An extensive linear array on the ground would contribute 


greatly to the study. 
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Figure Captions 
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1. Fluctuating values of total magnetic field obtained by a 
polar orbiting satellite after subtraction of a model field. 
(Cain and Hendricks, 1967). 

2. Magnetic effects of crustal anomalies observed at satellite 
altitudes over the United States (Zietz, e_t 1 970). 

3. Magnetic variations observed by a polar orbiting satellite 
as it passes over the equatorial electrojet (Cain and 
Sweeney, 1973). 

4. Anomalies in the equatorial electrojet signature illustrating 
day to day variability (Cain and Sweeney, 1973). 

5. Total field changes caused by satellite passage over the 
auroral electrojet (Langel and Cain, 1968). 

6. Transverse field variations observed by a vector field mag- 

netometer passing through a sheet of field aligned current 
flowing into the auroral oval (Armstrong, ^ 1974). 

7. High resolution ATS-1 magnetic field data showing effects of 
spacecraft field as discussed in text. Vertical scale 5 
gamma per division, horizontal scale 5.12 seconds per division. 

8. Power spectrum of the derivative of ATS-1 magnetic field data 
showing harmonics interference in transverse components at 
multiples of (1/5.12) HZ. 

9. Power spectrum of quiet ATS-1 field data showing harmonics of 
(1/81.92) HZ due to ramp in Z component. 
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Contour inttrvol 3 gonunai 

Residual magnetic map contoured from data taken by USSR satellite Cosmos 49. Contour interval is 5 7 . The altitude ranges from 261 to 4SS km. 
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Tho POGO data 



Examples of day and night traversals over western Pacific. Both day 
f^lid line) and night (dashed) show structure from 10® to 30®8 probably duo to 
magnetic anomalies. High *8houlders’ on each side of electrojet miniTnuTn in day 

trace is atypical. 



LATrrUOC (degrees) 

Magnetic imomaly around Bangui (Central African Republic) near position 
of deotrojet on both day and night passes. 
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Figure 3 
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Complex signature in field residuals during moderate magnetic disturbance. 



. Complex wiaiion with depression on one aide of jet exceeding that of 

the jet. 


Figure 4 
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GEOGRAPHIC LATITUDE 
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CL COMPARISON OF MODEL RESULTS WITH SURFACE DATA 


Disturtanwe fields from A polar current model compared with disturbances observed at OGO>2 and at surface observatories. 
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Tntroducti on 

A . Why should hydromaqneti c waves be measured in space ? 

HM waves are produced by plasma instabilities or currents 
far out in the magnetosphere and propagate to the earth's sur- 
face. For many hm waves there is reason to believe that propa- 
gation is f i el d-alti gned , or ducted. Thus, the waves seen at the 
earth's surface provide information about plasma processes going 
on at known points in space. The wave spectrum, polarization, 
and direction of propagation can be used to infer the nature of 
these processes. For example, the existence of Pc 1 waves 
(left elliptically polarized transverse waves of 5 sec period) 
has been used to infer the occurrence of the ion cyclotron in- 
stability of energetic protons in the magnetosphere. Furthermore 
the fact that one type of Pc 1 (pearl pulsations) shows disperion 
in its dynamic spectrum has been used to determine the L shell 
on which the waves are propagating, the equatorial density of cold 
plasma on this line, and the parallel energy of the protons re- 
sponsible for the instability. Another type of Pc 1, IPDP, is 
beginning to provide information about the point and times of 
proton injection during substorms as well as details about the 
drift of these particles. Other possibilities include the high 
latitude Pc 4,5 seen near the polar cusp, or Pc 3,4 seen near 
the plasmapause. These waves may eventually be used to provide 
information about the location of these boundaries as well as 
processes occurring on them. 

In summary, hm waves are diagnostic of plasma processes in 
the magnetosphere. The location at which these waves occur on 



the earth's surface reflects locations of important regions or 
boundaries in space. The properties of the waves reveal infor- 
mation about the plasma properties in these regions. The waves 
therefore provide the possibility of inexpensive monitoring of 
these magnetospheri c processes with ground observations. 

B. Current Feasibility of Monitoring Magnetospheri c Processes 
with hm Waves 

The general goal outlined above is not presently attainable 
for a numer of reasons. Mainly this is a consequence of our lack 
of understanding of the processes of wave generation and propa- 
gation to the earth's surface. This is, in part, the major justi- 
fication for an hm wave sensor on Space Shuttle. We presently 
have almost no understanding of the effects of the ionosphere on 
the transfer of waves. For some waves the ionosphere may act as 
a resonant cavity, increasing the amplitude of the incident waves 
much above their initial amplitude. Alternatively, the ionosphere 
may duct energy away from the point of entry. Almost certainly 
the polarization of the incident waves is greatly altered by the 
ionosphere. ^ 

These facts make it impossible to use current ground obser- 
vations to determine the point of entry of hm waves from the 
magnetosphere, or their polarization as it exists within the 
magnetosphere. Without this information, no clear inferences 
can be made about the locations or types of magnetospheri c genera- 
tion mechani sms . 

A major goal of an hm wave sensor on Space Shuttle would be 
,to determine the effects of the ionosphere on the transmission of 
hm waves to the ground. 


